ITEP-TH-22/12 



From elongated spanning trees to vicious random walks 

CN \ A. Gorsky 1 , S. Nechaev 2,3 , V.S. Poghosyan 4 , and V. B. Priezzhev 5 

1 ITEP, B. Cheryomushkinskaya 25, 117218 Moscow, Russia 

£h ■ 2 LPTMS, Universite Paris Sud, 91405 Orsay Cedex, France 

2 . 

3 P.N. Lebedev Physical Institute of the Russian Academy of Sciences, 119991 Moscow, Russia 

IT) 

^Institute for Informatics and Automation Problems NAS of Armenia, 375044 Yerevan, Armenia 

' I '' 5 Bogolubov Laboratory of Theoretical Physics, 

O 

D ■ Joint Institute for Nuclear Research, 141980 Dubna, Russia 

B 

Given a spanning forest on a large square lattice, we consider by combinatorial 
methods a correlation function of k paths (k is odd) along branches of trees or, 
equivalently, k loop-erased random walks. Starting and ending points of the paths 

s . 

are grouped in a fashion a k— leg watermelon. For large distance r between groups 

g 
o 

o 



of starting and ending points, the ratio of the number of watermelon configurations 
to the total number of spanning trees behaves as r~ v \ogr with u = (k 2 — l)/2. 



Considering the spanning forest stretched along the meridian of this watermelon, we 
see that the two-dimensional A;— leg loop-erased watermelon exponent v is converting 

""xT \ 

' into the scaling exponent for the reunion probability (at a given point) of k (1+1)- 

dimensional vicious walkers, v = k 2 /2. Also, we express the conjectures about the 
possible relation to integrable systems. 
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I. INTRODUCTION 



Enumeration of spanning trees on (i-dimensional grids is a classical problem of the combi- 
natorial graph theory According to the Kirchhoff theorem, the number of spanning tree 
subgraphs on a lattice is given by minors of the discrete Laplacian matrix A of this lattice. 
Considering a modified matrix A', which differs from A by a finite number of elements, one 
can define local correlation functions of spanning trees in any dimension. Besides of the 
combinatorial methods, in the two-dimensional case the conformal invariance can be used, 
which allows to evaluate the nonlocal spanning tree correlations of special topology. For 



From elongated spanning trees to vicious random walks 3 



instance, using a mapping of the g-state Potts model to the Coulomb gas, and taking the 
limit q — > 0, Saleur and Duplantier 2] have found the correlation function W k (r) of a fc-path 
"connectivity" in the spanning tree. The /c-path connectivity means that, given two lattice 
points A and B separated by the distance r (r ^> 1), one can find k paths along branches 
of the spanning tree from a vicinity of the point A to that of the point B. Such a topology 
is called in physical literature the fc-leg watermelon and each path along the branches is a 
loop-erased random walk (LERW). Schematically, the fc-leg watermelon is depicted in the 
figure Fig. [Th. In this work we will be mostly interested in the asymptotic behavior of W k (r) 
for large separations r>l. 

In the frameworks of the conformal field theory, the correlation function of k connected 
clusters tied together at the extremities separated by the distance r, has the following scaling 
behavior: 

W k (r) ~ r~ 2xk (1) 

where x k is the conformal dimension of the scalar primary fields, and k is the number of 
clusters. By mapping to the g-state Potts model with arbitrary q, one allows to see that 
each cluster of the Potts model in the limit q — > contains a path from one extremity to 
another. The values of scaling dimensions, x PtQ , have been computed as 

2h M , h m = [{m \ l)P f - m ^~ l (2) 



x P , q 4m(m + l) 

where h p>q are the weights belonging to the Kac table and p, q are integers in the minimal 
block, 1 < p < m — 1. The parameter m in ([2]) is fixed by the central charge c, 

c = 1 A r (3) 

m[m + 1) 

For c = —2 (i.e., m = 1) one has 

hM = (?ez|>!z! (4) 

and for (p,q) = (k/2,0) the large-distance scaling behavior of the A;-path correlation func- 
tion is 

W k (r) ~ r - 4hk ^° = (5) 

A combinatorial derivation of the exponent v = (k 2 — l)/2 for the fc-leg watermelon 
remains still an open problem. However, a "minor" modification of the problem makes it 
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solvable. Instead of the k— leg watermelon embedded into the single-component spanning 
tree, one may consider the two-component spanning tree and put the watermelon into one 
of these components. Then, for sufficiently large components, the partition function of k- 
leg watermelon configurations for odd k is given by the determinant of k x k matrix, whose 
entries are the Green functions A -1 of the discrete Laplacian A. 

The partition function of the 3-leg watermelon on the two-component spanning tree has 
been computed for the first time in [3J where the 3-leg watermelon, called the #-graph, has 
Deen used for exact determination of the height probabilities in the Abelian sandpile model 
4j, |5]. The k-leg generalization of the #-graph was the subject of the work [l0|, where the 
authors have derived a determinant expression for the fc-leg watermelon also embedded into 
one component of the two-component spanning tree. Calculating Wk(r) for few small values 
of k {k = 1, 3, 5), and guessing the structure of the series expansion of Wk(r) at r 3> 1, they 
have suggested the generic expression 



Wk(r) ~ r 2 Inr (6) 
where the logarithmic factor comes from the two-component topology of the problem (see 



e.g 



An exact derivation of ([6]) needs more careful analysis of the resulting determinant. 
In our work we show how the asymptotic correlation function ([6]) can be derived exactly 
from the corresponding determinant expression. 

The determinant expression of the fc-leg watermelon in the two-component spanning tree 
can be generalized to the situation when the spanning trees are elongated in some direction. 
Ascribing additional weights to the lattice bonds along some direction, one can obtain an 
asymmetric Laplacian A^ £ ^ which depends on the degree of elongation e (0 < e < 1/4), 
where e = corresponds to the isotropic lattice and e = 1/4 - to the fully elongated case. 
The Kirchhoff theorem, being applied to A^ £ \ gives the partition function of anisotropic 
spanning trees. 

The anisotropy affects the statistical properties of branches of the trees and, therefore, 
the loop-erased random walks constructed as paths on the elongated branches. At extreme 
stretching e = 1/4 along a diagonal of the square lattice, the loop-erased random walk 
becomes obviously the simple Bernoulli process with equal probabilities to take one of two 
possible directions. We consider the fc-leg watermelon statistics for e between and 1/4. In 
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the limiting case e = 1/4, one may expect that the loop-erased random walks constituting 
the watermelon behave like vicious walkers [6| and the watermelon geometry corresponds to 
the reunion at a given point in the space of k walkers after equal number of steps made by 
all walkers. The elongated k— leg watermelon is drawn in the figure Fig. [Tb. In contrast to 
non-stretched watermelon, it is natural to call it "zucchini". According to [6] and [7], the 
correlation function Wk(r) at r 3> 1 is 



Figure 1: a) Symmetric k-leg watermelon in 2D made of loop-erased random walks; b) Stretched 
watermelon ("zucchini") of k directed (1+1)D vicious walkers. 

If the behavior of the watermelon at extreme elongation is rather clear, the asymptotics of 
Wk(r) at small elongations is less obvious. The presence of the bias, produced by an external 
field, forces the loop-erased walk to follow the direction of the field, thus breaking the 
conformal invariance and making the two-dimensional system effectively (1+1) -dimensional. 
At small e, the external field is a perturbation of the logarithmic conformal field theory with 
c = —2, and determination of the exponent in (J7J) by the CFT methods becomes a non- 
trivial problem. So, the main goal of our work consists in examination how the watermelon 
exponent is converting into the scaling exponent for the reunion probability of k vicious 
walkers in the interval < s < 1/4. 

Another intriguing application of our work concerns the "interpolation between CFT and 
RMT". It is known that the isotropic LERW is described by the 2D SLE k with k = 2 
while the directed mutually avoiding ( "vicious" ) walks deal with statistics described by the 



W k (r) 



(7) 




x 
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random matrix ensembles. Changing the degree of extension, e, we pass from the c = —2 
CFT to the RMT ensemble. 



II. WATERMELONS OF LOOP ERASED RANDOM WALKS (LERW), 
SPANNING TREES AND WILSON ALGORITHM 

In order to make the contents of the paper as self-contained as possible, we outline in 
this section some connections between spanning trees and the loop-erased random walks. 

Consider a finite square lattice £ with vertex set V and edge set E. Select a tagged subset 
W (W C V), which is called "the set of dissipative vertices". Let V = [uq, Ui,u 2 , ■ ■ ■ ,u n ] be 
a collection of vertices (a trajectory) passed by a n-step random walk on C, started from 
Mo and ended upon hitting u n G W . The loop-erasure LE(V) = [70,71,72, • • • ,7m] of V is 
defined by removing loops from V in a chronological order. Note that the order in which 
we remove loops is important, it is uniquely fixed by the condition that the loop is removed 
immediately as it is created when we follow the trajectory V. The ordered collection of 
vertices LE(V) is called "the loop-erased random walk" (LERW) on C with dissipation W. 



Wilson {9] has proposed an algorithm to generate the ensemble of uniformly distributed 
spanning trees by LERWs. It turns out to be useful not only as a simulation tool, but 
also for theoretical analysis. The algorithm runs as follows. Pick up an arbitrary ordering 
V \ W = {vi, . . . , vn} for the vertices in C. Inductively, for i — 1, 2, . . . , N define a graph Si 
(5*o = W) which is a union of S^i and a (conditionally independent) LERW from Vi upon 
hitting Si-i. If Vi G SV-i, then Si = SV-i- Regardless of the chosen order of the vertices, 
Sn is a sample of uniformly distributed random spanning forests on C with the set of roots 
W. A spanning forest with one root is a spanning tree. A spanning forest with a fixed set 
of roots can be considered as a spanning tree, if we add an auxiliary vertex and join it to 
all the roots. When the size of the lattice tends to infinity, the boundary effects vanish, so 
we can neglect the details of the boundary. We do not distinguish between spanning forests 
and spanning trees in this case, assuming all boundary vertices are connected to an auxiliary 
vertex. 

Consider the set of k bulk vertices Ik = ■ ■ ■ , i>k}- The Wilson algorithm on £ 

with the dissipation W = DC U Ik generates k + 1-component spanning forest with roots 
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ii, %2i ■ ■ ■ i ik and an auxiliary vertex connected to the boundary vertices. Given k, and 
introducing another set of k vertices Jk = {ji, J2, • • • , jk}, we can enumerate all spanning 
forests, for which the vertices in Jk belong to the different components rooted at vertices in 
4- 

Construct now an auxiliary Laplacian on £ using the bridge trick (|3j). Specifically, define 
the unperturbed Laplacian of the square lattice, A, and the perturbed Laplacian, A', which 
differs from A by "defects" (bridges) with the weights — rj, 

4 if % 



4 if % — j 

— 1 if i,j are nearest neighbors i A- 
otherwise 



< 



J 

— 1 if i,j are nearest neighbors 
-77 if i,j is a pair (i s , j s ), s = 1, 
otherwise 



One can check that the combination 



lim — detA' 



,k 
(8) 

(9) 

rj— >oo Tj 

gives the sum over k + 1-component spanning forests T having k bulk roots in the lattice 
points Ik = {ix,iz, . . . , ik}- Each forest T enters into the sum with the sign (— l) c( - T \ where 
c(T) is the number of cycles which appear in T due to adding bridges (i s ,j s ), s — 1, 2, . . . , k 

- see . 

Following Ivashkevich and Hu jlOj, we chose sets Ik and Jk as zigzags (or "fences") with 
odd k - see Fig. [2^,. Then, all possible configurations of cycles on the bridges contain 
either 1 or k cycles - see Fig. [3] and, therefore, (— l) c(T ) = —1. For fixed k and the large 
distance between sets Ik and J/-, the bunch of lattice paths from Jk to Ik has a form of a 
watermelon. Connecting the neighboring points in the set Ik (and in the set Jk), we obtain 
a two-component spanning tree with the watermelon embedded into one component. Thus, 
we have 



lim — det A' = number of fc-leg watermelons 



(10) 



rj— >oo fj n 

The matrix A' can be represented as A' = A + B, where the matrix B has nonzero 
elements —n for pairs (i s ,j s ), s — 1,2, . . . ,k. To evaluate the determinant A', we use the 
known formula jl2| 



det A / 
detA 



det(I + BG), 



(11) 
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Figure 2: The k nonintersecting loop-erased random walks begin in the set Ik and end in the set 
Jfc: (a) Sets Ik and (k is odd) chosen as parallel zigzags on the square lattice; (b) Sequential 
zigzags Ik and Jk oriented along the horizontal axis. 





Figure 3: Possible closures of zigzags by loop-erased random walks. 

where I is the unit matrix, G is the Green function G = A^ 1 and the nonzero part of the 
matrix B is the k x fc-matrix. Inserting (TTTj) into (jTUl) . we obtain the ratio between the 
number of two-component spanning trees containing the fc-leg watermelon and the total 
number of spanning trees on lattice C: 



W(h, J k ) =detG(4, J k 



(12) 



where rows and columns of matrix G(Ik,Jk) are labelled by indices h,i2, ■ ■ ■ ,h and 

31J2, ■ ■ ■ ,jk- 
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III. ISOTROPIC LERW WATERMELON 



The correlation function W(Ik, Jk) at large separation, r, between the sets Ik and J k 
(k is odd) takes an asymptotic form Wk(r) ~ r _I Tnr [h]], where v is a universal lattice- 
independent exponent. To find we have to specify the sets Ik and Jk on the square lattice 
and consider then the continuous limit r — > oo. Following (|T2l) we should define the matrix 
G(Ik, Jk)- 

The Green function G x<y on a square lattice is the solution of the equation 

+ G X,y+1 ) + G X: y — 5 X fl5yfi (13) 

The integral form of G x>y is 

Gl „ = ^ T *. [' n, ^)^y) (14) 

8n 2 J_ 7T J_ 7T 2 - (cos a + cos/3) 

For a specific configuration of initial and final points, we can fix: x = x(ji) — x(ii), y = 
y{ji) —y{ii), where the pairs (x(ik), y(ik)), ( x (jk),y(jk)) label the coordinates of the vertices 
ikiik i n the sets Ik and Jk respectively. Define r = a/x 2 + y 2 - At r = the Green function 
G XjV is singular, and at large distances, r>l, the function G X)V = G(r) has the asymptotic 
form 

G(r) = G(0 )-^l„r-£-?£? + ... (15) 

where 7 is the Euler constant. 

Taking the sets of initial (Ik) and final (Jk) points as zigzags depicted in the Fig. [2^ 
and fixing the coordinates of left ends, i\ e Ik and ji G Jk as x(ii) = 0,y(i\) = and 
x (ji) — x i y(ji) — we can easily define the coordinates of the vertices m in zigzags Ik and 
Jk- So, we have 



m ,. s m — l 

(16) 



for m even, and 



, . , m — 2 m 

x(] m ) =x+ , y(j m ) = y - — 



m — 1 N m — 1 

x ( ? m) = — - — , y(*f 



2 2 (17) 

m — I m — 1 

zUn) = aH — ^ — , y(j m ) = y o — 



for m odd. 
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Now we can explicitly write the expression for the determinant W(Ik, Jk) — det Jk) 
in f lT2|) for specific position of zigzags Ik and Jk (k is odd) shown in the Fig. [2j 

G \ 



W(I k , J k ) = det 



'x,y 



G x ,y~l G x+ i : y-i G x+ i t y-2 ' ' ' G x +k y y-k 



G x —\^y (*x—l,y—l G Xy y—\ G Xy 



y-2 



Gx— G x — 1,2/ G X) y G X y—\ 



Gx-2,y+l G 



,y+l ^x-2,y 



y -'x—\,y ^x— 1 



(18) 



The generic form of the entry in the position (m',m) of the matrix in (fTBl is G XmmhVmm ,, 
where x m , m / = x(j m ) - x(z m /), y m , m / = y(j TO ) - y(i m >), {m, m'} = l,...,k. 

The explicit expression f|T8|) is useful for the determination of the exponent v for small 
k = 1,3,5, .. .. However, the origin of the general expression v = (k 2 — l)/2 remains still 
hidden. To reveal it, we note that the bridge construction used in the previous section in 
the watermelon definition does not depend on the relative orientation of the sets Ik and Jk- 
Then, we can chose a more convenient orientation shown in Figj2b. To make the round 
along bridges possible, in the latter case we should insert bridges (i s ,jk+i- s ), s — 1, 2, . . . , k 
instead of former ones (i s ,j s ), s = 1, 2, . . . , k. 

The asymptotic analysis of Eq. fll2l) does not need exact locations of the points ii, . . . , ik 
and ji, . . . , jk in the sites of the square lattice. Instead, we can chose for i lf . . . , i k and 
ji, . . . ,jk two sets of horizontal coordinates < U\ < u 2 < ■ ■ ■ , Uk < const and r + v x < 
r + v 2 < ■ ■ ■ < r + v k with < v , < const, i = 1, . . . , k. The vertical coordinates of points 
ix, . . . , ik and ji, . . . ,jk will be for odd indices and £ > for even ones. Given coordinates 
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{ui} and {v^, the r-dependent part of the determinant in (fl2|) has the asymptotic form 
ln[(r - Mi + vx) 2 } \n[£ 2 + (r - u 2 + vx) 2 } \n[(r - u 3 + Vi) 2 } ■■■ ln[(r - w fc + v i) 2 ] 

ln[£ 2 + (r - Ul + w 2 ) 2 ] ln[(r - u 2 + v 2 ) 2 } \n[(£ 2 + (r - u 2 + v 2 ) 2 } ■■■ 
ln[(r - m + v 3 ) 2 ] ln[£ 2 + (r - u 2 + w 3 ) 2 ] hi[(r - u 2 + f 2 ) 2 ] 

ln[(r - m + w fc ) 2 ] ln[(r - u k + w fc ) 2 ] 

(19) 

The matrix elements in Eq. (fl9|) are of two types: ln[(r — Ui + Vj) 2 ] for odd-odd and even- 
even indices i and j and ln[£ 2 + (r — ttj + u,) 2 ] for even-odd and odd-even indices. For large 
distance r, we have r ^> u^, r ^> v i: i = 1, . . . , k and r ^> £. We rewrite the matrix elements 
as 21nr + 21n[l — Uj/r + fj/r] and 21nr + ln[(£/r) 2 + (1 — u-Jr + Vj/r) 2 ] and expand in 
powers of £/r, Uj/r and fj/r. 

Consider first the case £ = 0. The leading term of the expansion of the determinant ffl9|) 
is a totally antisymmetric polynomial of the form 

k k 

Hiui-uj) Hivi-Vj) (20) 

j>i j>i 

Each Ui and brings the factor 1/r, so that the leading term of the expansion is of order 

Consider now the case £ > 0. The sets and </& split into two subsets each, having vertical 
coordinates £ = and £ > correspondingly To distinguish between different subsets, we 
supply coordinates Ui and Vi for even % with hats: Ui and and re-enumerate each subset in 
increasing order. We obtain four sets: I m = {ui, u 2 , . . . , u ni }, I n2 = {ui, u 2 , ■ ■ ■ , u n2 },J ni = 
{vx,v 2 , • • • and J„ 2 = {vi,v 2} . . . ,< 2 }, where m = (A; + l)/2 and n 2 = (k - l)/2. 

The resulting determinant is a totally antisymmetric function with respect to arguments 
Ui, u 2 , . . . , u ni , Mi, u 2 , . . . , i4 2 , t>i, f 2 , • • • , u»n an d ui, v 2 , . . . , Vn 2 separately. Therefore, the 
first survived term in the expansion of the determinant has the form 

ni ni ni m 

Y[( Ui - Uj) Y[(Vi - Vj) ]J(Ui ~ Uj) Y[(Vi ~ Vj) (21) 

j>i j>i j>i j>i 
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and the contribution from u, u, v, v to the leading term of the expansion is of order r In r, 
where 

k 2 1 

^ = n 1 (n 1 -l) + n 2 (n 2 -l) = — -k + - (22) 

Besides the horizontal coordinates u, u, v, v, we have to take into account the vertical coor- 
dinate £. The minimal order of £ in the expansion of the determinant is 2n 2 — k — 1. Indeed, 
assume that the order of £ is less than 2n 2 . It means that at least one row or column of the 
matrix f|T9|) is free of £, and as consequence, at least one element of the set J n2 or J n2 should 



be transferred to I ni or J ni . Then, we obtain instead of (122]) 

, 774(774 - 1) n 2 (n 2 - 1) »i(wi + 1) (ra 2 -2)(ra 2 - 1) A; 2 ,5 
^= 2 + 2 + 2 + 2 = y- A;+ 2 >/i (23) 

Therefore, the minimal order of £ is k — 1 and the resulting exponent in the expansion of 
the determinant is 

k 2 - 1 

v = n + k - 1 = — - — (24) 

The sets 7^ and Jk are particular cases of the more general sets Ik = I ni U ^2 an d Jk = 
Jm U ^"25 so the exponent in (JB]) coincides with v. 



IV. ELONGATED LERW WATERMELON 

A. Massive Green functions 

The determinant expression (jT2l) for the correlation function of a bunch of k loop-erased 
random walks remains valid even for stretched spanning trees. When the external elongating 
field is applied, the theory becomes massive. That induces an effective lengthscale in the 
ensemble of biased loop-erased random walks. 



Let G Xj y(e) be the lattice Green function of a two-dimensional random walk, with the 
drift controlled by e. The function G XjV (s) satisfies the following equation (compare to ( TT3|) ) 
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The substitution (|26|) allows to rewrite ([25]) in the symmetric form: 

- a$Q x -i t y(a) + Q x>y -\(a) + Q x+hy (a) + Q x>y+ i{a)\ + Q x , y (a) = A~ {x+y) l 5 xfi 5 yy0 (27) 



where a = a/1 — 16e 2 . By the Fourier transform we get (compare to ([T 

^ , s 1 /' 7r , cos(ax) cosiBy) 

Q*,„(a = 5-a / «W ^ o 7 m 28 

8n 2 J_ w J_ w 2 — a(cosa + cosp) 

Thus, the Green function G X)y (e) with the drift e is: 

1 fl+ie\^ M P P cos(ax) cos(/fy) 



G x v (e) = — - — ' / da d/3 v ' ^ (29) 

M ' 8tt 2 Vl - 4^7 J_ n P 2-v / T^W(cos« + cos/3) 1 ' 

Rewrite G x>y (e) as 

Gf S)tf (e) = \ (\^J 2 dre- 2T I x (yft - 16e 2 r) /„ (V 1 - 16 £ 2 r) (30) 

where I x (...) and I y {---) are the modified Bessel functions of integer orders x and y re- 
spectively It is instructive to re-express the equation ( 1301) in terms of the hypergeometric 
function 4F3: 

fl+x+y 2 + x + y 1+x + y 2 + x + y 2 X 
4F3 , , , ; 1 + x + y, 1 + y, 1 + x; 1 - 16e 



The equation ( 1311) . being an exact expression for the Green function of the random walk on 
in open square lattice (x, y) with the drift along the diagonal x = y, allows one to analyze 
the asymptotic behavior of G Xjy (e) (at r = a/x 2 + y 2 ^> 1) for e — > \ (extreme elongation) 
and for e — > (no elongation). 

B. Extreme elongation (s — > |) 

To consider the strong elongation, we extract the asymptotic behavior of G Xty (e) at e — > j 
and plug this expression into the determinant formula ( 1T21) . The expression ( 13T]) allows to 
find the desired asymptotic behavior of G x>y (e) easily. Introduce the new variables: 

x + y x — y 

- 1 = ^r ; 2l = ^r (32) 
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and rewrite the Green function (13TT) in terms of z\\ and z±: 



4 V4 / r(i + 3r^i)r(i + 



' F3 (l + Tr 1 + ^'5 + 7=2 ' 1 + 71 ; 1 + ^" 1 + ZJ 1T ■ 1 + 1 " lfe2 ) 

(33) 



At e — >■ I we have 



G " » = rfiTSrn^L^ [1 + (1 " i£)Rl + (1 " 4C)2R2 + 0(1 " 4£,3] (34) 
where -Ri = Ri(z\\,z±) and i?2 = R2{z\\,z±) are rational functions of z\\ and 

itl^||,3jj - 2(2+2^,1+4-4) 

(35) 

6z|[+v / 22j|(44+35)+222(24+184+73)-2v / 2z||(4- 22 4- 64 )+ 8 (' z i+ 10 ) 
" ' ~ L ' 16(2jf+6v / 22|-22j|(4-13)-6v / 22|| (4- 4 )+-4- 10z i+ 16 ) 

Denoting a; = \ — e and applying the Stirling formula to lnG^ 2± at z\\ ^> 1, we get the 
following asymptotic behavior at large separations: 

1 + 2, 4 -y _ 1 1 1 + 44 + Ml - 44) + 2^(7 - 44) Q / 1 \ 3/2 
2 11 4^,^(1 + 20; + 6a; 2 ) W 

(36) 

Keeping only the first order terms in lo = \ — e, we get from ( 136]) the following asymptotic 
behavior for the Green function Gl D z± at large separations: 



•2 



- ^e"^' - ^ 7=F^ e~ <*o (37) 



^2(l+2u)z|| ^ c v^z 

2 7 /4 v /7f^ ~ 2^4^(1 _4w)7T2|| 

Thus, we see that (|37|) has a form of a Gaussian distribution for a one-dimensional random 
walk of length z\\ and the transversal coordinate z± for any strong elongations < u = 
2 — e <C 1. The separation distance z\\ between two ladders shown in the Fig. [2] can be 
interpreted length of a directed trajectory. 

For e = |, the Green function G^P Z± becomes the probability to visit point x = z± of 
the one-dimensional lattice {x G Z 1 } by the symmetric random walk after z\\ time steps, 
provided that it is located at the origin x = at the initial moment of time z\\ = 0. The 
sets Ik and in the determinant formula f lT2|) can be interpreted in terms of variables 
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zu and zj_ simply as groups of k neighboring sites of the one- dimensional lattice taken at 
the moments of time t — and t — z\\. Then, formula ffl2|) takes the form of the Karlin- 
McGregor formula (known also as the Lindstrom-Gessel-Viennot formula and used by Huse 
and Fisher for calculation of the reunion probability of the vicious walkers). The result 

n q 

obtained in [a, [7[ for large number of steps t is given by (J7J) where r — zn — t. 



C. Arbitrary elongation e > 



The direct expansion of fl3Tl) for small e and large x, y meets essential technical difficulties. 
Taking into account that for i>1 and y S> 1 the main contribution to the integral in (j29|) 
is given by values \a\ < | and \(3\ < K Thus, at large separations r> 1 and at small e, we 
can use the saddle-point approximation. To proceed coherently let us replace the discrete 
recursion relation ( 12"5"|) by the continuous one 

- 2e (d x + d y ) G x>y {e) + \ (d 2 xx + d 2 J G XjV (e) = 5 x>y (38) 

Note that the equation (|38|) is valid for any r and e. The solution to ( )38|) is 

e &£(x+y) [-co roc e i(xa+y/3) 



-—L da L i(i *> + H* + p) (39) 

The expression f )39|) can be easily evaluated exactly. Introducing the polar coordinates 

a = p cos (ft; (3 = p sin 

rewrite (13"9"|) as 

g4e(a;+y) ^00 i-2tt ^ip{x cos cf>+y sin cf>) 

Since 

x cos + y sin = \/ x 2 + y 2 cos(0 — 7) 
where 7 = arccos -, x , we get: 

yx 2 +y 2 

p ie{x+y) poo J Q (p^/x 2 ~Ty 2 ) 0p^ x+ '^ / r- , \ 

G x , y {e) = / pdp ) a , l2 J = K 4v / 2e v / ^+7) (41) 

Rotating the coordinates by 7r/4, i.e. passing to (^|| , ) defined in fl32|) . we arrive at the 
final expression 

G x>y (e) = —e^Ko {a^e^J ' z 2 + z\) (42) 
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To extract the asymptotic behavior of G Xty (e) at any e > and large Z\\, we should first 
expand logarithm of (T42]) ar z\\ ^> 1 and then consider the elongations with any e > 0: 



In G x , y (e) = - In (2 3 /^) - I In *„ - + O ^ ] (43) 

Thus, we get from ( 1431) at zn ^> z± ^> 



4ez 



2 



In Appendix we demonstrate that similar expression can be obtained for stretching along 
the x-axis. To demonstrate the generality of the result, we perform the computations for 
the discrete model. 



V. ELONGATED LERW WATERMELON VS VICIOUS WALKERS 

The explicit expression for the determinant (1121) is 

W(I k , J k ) = S(n)Gi^ 1 )G 2t7T (2) ■ ■ ■ G kj7T ( k ) (45) 

where 5(ir) is the sign of permutation 7r and G m ^ n is the Green function G x jV defined 
in (USD. 



The expressions (1371) and (1441) allow us to write (1451) in the form 



W(I k , J k ) = £ (46) 

The diffusion constant D = D(e) depends on elongation e and defines the mean-square 
displacement ((zj — -2°) 2 ) = Dr in the transverse direction for the "time" r. The coordinates 
Zi, . . . z k and zf, . . . z\ of vectors z and z° denote the positions of points z'x, . . . i k and ji, . . . j k 
measured along corresponding zigzags. The expression (146]) coincides with that used by Huse 
and Fisher (see Appendix of [3]) for the derivation of asymptotic reunion probability of k 
vicious walkers. 

The computation of (146]) is based on the analysis of antisymmetric properties of variables 
Zi, . . . z k and . . . z k as well as on the computation of ([19]) was based on the analysis of 
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antisymmetric properties of variables ui, ...u k and Vi, ...v k . It is instructive to compare 
procedures leading to isotropic and elongated asymptotic behavior. 

Using the invariance |vrz°| = |z|, one can write (1461) as 



| z | + 1 z 



0,2 



W(I k ,J k ) = U k ( Z y)^£ ¥J - 2 (47) 

where 

U k (z,z°)= (48) 

TT<=S k 

Since the length of zigzags L is much smaller than the separation between clusters I and J, 
one can expand the exponential in powers of Each term of the expansion is a homo- 

geneous polynomial in k variables or, equivalently, The lowest order of polynomials 
is of degree k(k — l)/2 due to the requirement of total antisymmetry. Each Zi or 2? brings 
the factor 1/r, so the asymptotic behavior of U k is U k ~ l/r fe2 / 2-& / 2 . Combining this factor 
with l/(2irDr) k ' 2 in ([47]), one gets the asymptotic reunion probability (J7J, Wfc(?") ~ r 2 . 

The first difference between the isotropic and elongated watermelons is in the depen- 
dence of the determinant expansion in m and Vi for isotropic watermelons and in Zi and 
for stretched ones. The arguments and zf enter into the expansion pairwise, while the 
expansion of f|T9|) depends on U{ and Vi independently, so both and Vi bring the factor 1/r 
each. 



The second difference is that matrix elements of ( 1451) are uniform, while those of ( fl9l) 
contain the additional argument £ depending on the parity of indices. Each £ increases the 
power of 1/r in the expansion by one (see the explanations at the end of the Section III). 

Finally, the determinant f|46l) contains a "diffusive" factor 1/ 1 (2nDr) k l 2 originating from 
the asymptotic form of the Green function (1371) . Being collected all together, these factors 
lead to different laws (jSJ) and (J7J for isotropic and elongated watermelons. 

VI. DISCUSSION 
A. Reunion exponents for LERW and VW watermelons 



Using the combinatorial approach based on exact enumeration of spanning trees on a 
square lattice, we have reproduced the Karlin-McGregor (or Lindstrom-Gessel-Viennot for- 
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inula for vicious random walks, or Slater determinant for free fermions) from the determinant 
expression for reunion probability of loop-erased stretched watermelons. 



Let us recall the Karlin-McGregor formula [13[. Consider k (1+1) -dimensional Markov 
chains Ci(*)> •••> Ck{t) with Q(t = 0) = x^ and Q(t = T) = xf\ for all i = 1, k. Assume 
that the time, t, is homogenous and any two walks cannot touch each other. Then the 
probability P T {d(t) < £ 2 (t) < ••• < (k(t)}, (0 < t < T) to have a bunch of k uncrossing 
random walks each of length T with starting points {C(0)} and ending points {((T)}, is 
given by the following expression (Karlin-McGregor formula) 

MCiit) < Ca(*) < - < Ck(t)} = det [P T (0(0), CjCO)i<jj< n ] ( 49 ) 

where Pt (C«(0), Cj(T)) is the probability for a single (l+l)-dimensional random walk of 
length T to have two extremities at the points d(0) and Q(T) (1 <i,j < k). 

In the original determinant expression (fl2l) (for e = 0) we have no any time scale in 
the system since the building block G X)V of the matrix W in ( fl2i) is the logarithmically 
decaying propagator of the 2D free field. However, when the bias e > is switched on, the 
theory becomes massive, what signals the occurrence of a finite length scale in the system, 
i ~ e~ l . The corresponding partition function, G XtV , becomes Gaussian at strong elongations 
controlled by the parameter e (see Eq. (l37|) ). Thus, the equation (fl8|) initially written for 
a A;-leg symmetric watermelon of two-dimensional loop erased random walks, converts at 
(0 < I — £ 1) into the Karlin-McGregor formula for a bunch of k (l+l)-dimensional 
vicious random walks. 

The value of the reunion exponents v = for the 2D fc-leg watermelon transforms 
at strong elongations into the (1+1)D value v — Note, that v differs by the constant 
— | from a reunion exponent v' = of the (1+1)D vicious random walks. The difference 
between v and v' deals with the fact that the value v obtained from Karlin-McGregor 
formula imply that the end points of the bunch meet in some specific point, while to get v 1 
we should allow for a reunion in every point in the space, i.e. intergrate over the position of 
the end point of the bunch. 
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B. Geometry of elongated watermelons and dimers 



The elongation of spanning trees and the passage from isotropic to directed watermelons 
made of loop-erased random walks, has a transparent geometric interpretation and can be 
understood in terms of configurations of dimers on square and honeycomb lattices. 

It is known (see, for example, [ijj]) that there is a bijection between the configuration of 
dimers on a square lattice Z 2 and the a spanning forrest on, say, an odd-odd sublattice Z^ 
of Z 2 . This bijection is set by the construction described in [14J. In brief, this construction 
is as follows. Consider the dimers located on the odd-odd sublattice Z^ only. Replace the 
dimers by vectors whose starting points are in the vertices of this sublattice. Continue these 
vectors until meeting the next neighboring dimer (vector) - as it is schematically shown in 
the Fig. HI Erasing arrows, we get on Z^ the spanning forrest, which is in bijection with 
the initial configuration of dimers. The configuration of dimers on the even-even sublattice 
can be uniquely restored since dimers fully cover the square lattice Z 2 without common 
points and empty vertices. 




II 



it 



Figure 4: Configuration of dimers on the square lattice and corresponding configuration of sym- 
metric spanning tree on the odd-odd square sublattice. 

The same algorithm we can apply to dimers on a honeycomb lattice (drawn in a form 
of a brick lattice). This establishes the bijection between dimers on a brick lattice and 
elongated spanning forest on an odd-odd square sublattice Z^ - see Fig. |5j Thus, the fully 
elongated spanning forrest (along the vertical axis in Fig. |5]) is in one-to-one correspondence 
with configurations of dimers on the honeycomb lattice. According to the correspondence, 
discovered by [15] (we have reproduced it schematically in the Fig. EJ), the configurations 
of dimers on a honeycomb (brick) lattice is in bijection with ensemble of directed mutually 
avoiding lattice random walks ("vicious" walks). 
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Figure 5: Illustration of bijection between dimers on the honeycomb lattice and fully elongated 
spanning forrest along y axis on an odd-odd square sublattice. 

As we have seen in the previous Section, the extreme elongation of loop-erased watermel- 
ons on the square lattice brings us to statistics of vicious walks extremely stretched along 
the diagonal of the square lattice. The geometric connection between spanning forest on a 
square and on a honeycomb lattices is consistent with this results (up to the rotation of the 
direction of extension by 7r/4). 



C. From CFT to RMT 



We would like to pay attention to the following observation. We see that the critical 
exponent v^n for e = defines the scaling behavior of the fc-leg correlation function of 
simplectic fermions described by c = —2 logarithmic conformal field theory and the reunion 
probability is given by the determinant expression (|T8l) . For < e < \ the theory is not 
conformal anymore because of the presence of the mass. At strong drifts, e ~ \ the system 
becomes (l+l)-dimensional with the critical exponent v — ^ and the reunion probability 
is described by the Karlin-McGregor formula which can be straightforwardly derived from 
( fT8l) by taking e — >■ \ limit. 

The interpolation between the determinant (j!8p and the Karlin-McGregory ( 149]) formulas 
as e is changing from to 4, allows us suggest the possibility to pass in the framework of 
unified description from the logarithmic conformal field theory to the random matrix theory 
since at extreme elongations the bunch of directed watermelons (vicious random walks) 

n 

represents the bunch of interacting Dyson random walks [16( described by the Gaussian 
Random Matrix Ensemble. 
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D. From LERW to VW via integrability 



According to [17], the isotropic LERW in the continuum limit is described by SLE^ 



corresponding to the conformal field theory with the central charge c = 1 



21 



For k = 2 one arrives at logarithmic CFT with c = —2. Let us use now the relation between 



SLE2 and the trigonometric Calogero system found by Cardy in 171 ] . He has shown that 
the stochastic radial SLE evolution is governed by the trigonometric Calogero Hamiltonian 

lpS 2 P(P-2) . _ 2 Xj - x k 

^ = -5L^ + nr- z. sm —5— ( 5 °) 



2^<9^ 



where — |. For k = 2 the coupling constant corresponds to the symplectic ensemble, 
/3 = 4. The n-particle Calogero wave function depends on the position of the SLE roots on 
the circle and corresponds to the probability in the LERW model. Simultaneously, the wave 
function can be identified with the particular conformal block in the c = —2 model 



¥col(*L, ■■■Zn) 0C< $(0)y 2 ,l(^l) • • • V 2 ,l{z n ) > 



(51) 



where V^i is associated with the states degenerate at the second level. 

The observation by Cardy could be used for the interpolation between different stochastic 
systems. Indeed, there exists an expected correspondence between the stochastic processes 
and various quantum Hamiltonian systems, so the Calogero model is not an exception. 
Since no elliptic parameters (or, equivalently, nonperturbative effects) involved, the most 
general relevant integrable many-body system is the trigonometric Ruijsenaars model which 
is the relativistic/discrete generalization of the Calogero system. The wave function of 
the trigonometric Ruijsenaars model is identified with the Macdonald polynomial. The 
corresponding stochastic counterpart, the Macdonald process, parameterized by (q,t), has 
been introduced recently in 20]. 

There are several degenerations of the Macdonald process 20|. The simplest one q — > 1 



brings us back to the trigonometric Calogero model. A more involved, so-called Inozemtsev 
scaling limit 2l|, provides the relativistic Toda system or, the g-Whittacker process, in 
the terminology of f2o|. Further deg eneration yields the usual Toda system or Whittacker 



process 



also 



25|, 



introduced by O'Connell 22|, [23), which has been found in several models other (see 



261]) . Let us emphasize that the coupling constant t in the rational Ruijsenaars 
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model fixes in general situation the central charge of the Virasoro algebra similar to the 
Cardy's case. The list of mentioned degenerations is schematically depicted in the figure 
Fig. H 
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V J 
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on a cylinder 
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Figure 6: The flowchart of possible degenerations of the trigonometric Ruijsenaars model. 



Turning back to the anisotropic lattice LERW model considered in the paper we could ask 
about its integrable counterpart. Since the isotropic continuum LERW correspond to the 
trigonometric Calogero model at fixed coupling constant, we should find its generalization 
in two direction because we have now anisotropy and the lattice size parameters. 

Conjecture. Since the lattice parameter implies the relativization of the Calogero model, 
the anisotropic LERW is mapped onto the relativistic Calogero model with an arbitrary 
value of coupling constant (i.e. of central charge). 

What could be the evidences in favor of such identification? 



• On the one hand, looking from the LERW side, we know that LERW model is closely 
related to the dimer model on the squire lattice. The integrable system associated to 
dimers, is the relativistic Toda chain |27| , which itself is the limit of the trigonometric 
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Ruijsenaars model, as we expect. However yet it is not clear how the Inozemtsev limit 
is set up in terms of dimers. 

On the other hand, we expect that the integrable counterpart of discrete vicious walks 
is the model of free relativistic fermions involved into the Chern-Simons action. Indeed 



it has been found in 



18 



19] that the reunion probability in VW model is identified with 
the 2D Yang-Mills theory on the sphere. It is well-known that lifting YM theory to the 
perturbed Chern-Simoms with the same group gives the relativization of the system, 
where the level of the Kac-Moody algebra provides the relativization/discretization 
parameter. The CS theory degenerates to 2D Yang-Mills in the k — > oo limit. 

It is useful to have in mind that the trigonometric Ruijsenaars model emerges from the 
perturbed Chern-Simons model on the torus with the marked point [29], while the trigono- 



metric Calogero model corresponds to perturbed 2D Yang-Mills theory on the cylinder [28]. 
Hence the stochastic process is formulated in terms of the particular observables in the per- 
turbed topological gauge theories. The coordinates of the particles involved are obtained 
from the holonomies of the gauge connections over the cycles upon the T-duality trans- 
formation. In particular, coordinates of the wave function of the Calogero model, which 
can be visualized via the cross-section of the watermelon in the middle, corresponds to the 
holonomies in the 2D Yang-Mills model. Note that we start with the discrete model in 2D 
hence at the extreme elongation we get the (1+1) system with the discrete time dynamics. 

The relation with the integrable models suggests another analogy which could be of some 
use. It is known that the Calogero model is related to the Quantum Hall Effect (QHE) at the 
continuum. The external magnetic field defines the coupling constant in the Calogero model. 
The type of Calogero model depends on the geometry of the QHE system. Suppose that 
we are in the limit of the strong magnetic field. Then the degrees of freedom are localized 
at the lowest Landau level and the 2D system becomes effectively (1+1) dimensional where 
the time is identified with the radial coordinate on the plane. This reduction is very similar 
to the one considered in the paper, however the decomposition of 2D into (1+1) is different. 

It would be interesting to see if two different conformal behaviors similar to discussed 
above could be found in QHE. Since the interpolation parameter is the magnetic field 
strength, we could ask what happened when we change it. The Calogero system describes 
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the gapless boundary modes of the QHE droplet, corresponding to the CFT with the central 
charge fixed by the filling fraction (or, equivalently by the magnetic field). What is the c=-2 
counterpart in the case of "magnetic elongation" ? Fortunately the answer is known: and the 
logariphmic conformal field theory describes the transition between the two QHE platoes. 
Hence, when the strength of the magnetif field is increasing, one evolves from one minimal 
CFT to another with the LCFT critical behavior. It seems very interesting to pursue further 
this analogy further. 
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Appendix A: Extension along an axis: Saddle point for discrete model 

Consider the Green function G of the Laplacian A on the anisotropic square lattice 
defined as a solution of the equation 



Here the parameters of anisotropy (drift) p±, q±, P2 and q<i are positive and satisfy the identity 
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^G Xt y — P\G X -l >y — 



QlG X+ l : y — P2G X: y-\ — 




(Al) 



Pi + qi + P2 + <?2 = 4 



(A2) 



It is easy to see that the integral representation of the Green function is 




(A3) 



Let us introduce the notations 




(A4) 
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Changing the variables a and f3 

e ia -> 7 ie iQ , e ip ->• 72 e i/3 (A5) 

we get 

G ^y = o 2 x y da dp - — - (A6) 

8vr 2 7f 7| ^ 2 - cos a - £ 2 cos /3 

Integrating over a, we arrive at the following expression: 

1 f ( 2 " ^ 2 cos ^ " V^-^cos/?) 2 -^ 2 ) * cos(y/3) 

G -» = 2vr(67i) X 7l io ^ ^-^cos/?) 2 -^ (A?) 

The asymptotics of the Green function for fixed finite y and x ^> |y| can be obtained 
using the saddle point approximation. To proceed, expand the exponent in (IA7I) near the 
maximum (3 = 0. 



2 - 6 cos $ - J (2 - ^ cos/3) 2 - £ 2 = 2 - 6 - V(2 - 6) 2 - £? 



X - P {- 2v /(2 6 ^- a } (l + 0( ^ <A8) 
The denominator of the integrand can be expanded in series in P and integrated term by 
term. For the main contribution, we have 

(2-6-v / (2-6) 2 -^i 2 ) X ' r - ^ ■ 

~ ^ 7 / d/3 e V< a -«>> a -«? cosiyP) (A9) 

^ 2tt (£ l7 i)%V(2 " 6) 2 - £ 

which gives 



G x , v ^ — its — ( — - — V — - — -I > x »m ( A1 °) 

Next terms of G x>y can be found by keeping higher powers in p and integrating sequentially. 

Assuming in (IA9I) that P2 = <?2 = 1, i-e. £2 = 72 = 1, and pi = 1 + 6, q% = 1 — S with 

5 > 0, we get 



G SiW ~ - — ^ ^ e v — , a; > (All) 



1 _j£5 

2v / 2vrv / 5 v / ^ 

Next terms of the asymptotic expansion of the Green function for the symmetric walk in 
vertical direction (p 2 = Q2 — 1) and fully elongatated in horizontal direction (5 = 1) are as 
follows 

G i_fr_£ + i>'-V + i_V-8V + W + i^ 7 
V 2x 8x 1921 / 
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